home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DSPSL.z / DSPSL
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSPPPPSSSSLLLL((((3333FFFF))))                                                            DDDDSSSSPPPPSSSSLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSPSL   - DSISL solves the double precision symmetric system A * X = B
  10.      using the factors computed by DSPFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE DSPSL(AP,N,KPVT,B)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAAPPPP DOUBLE PRECISION(N*(N+1)/2)
  20.         the output from DSPFA.
  21.  
  22.      NNNN INTEGER
  23.         the order of the matrix  A .
  24.  
  25.      KKKKPPPPVVVVTTTT INTEGER(N)
  26.         the pivot vector from DSPFA.
  27.  
  28.      BBBB DOUBLE PRECISION(N)
  29.         the right hand side vector.  On Return
  30.  
  31.      BBBB the solution vector  X . Error Condition
  32.  
  33.      AAAA division by zero may occur if  DSPCO  has set RCOND .EQ. 0.0 or  DSPFA
  34.      has set INFO .NE. 0  .  To compute  INVERSE(A) * C  where  C  is a matrix
  35.      with  P  columns
  36.         CALL DSPFA(AP,N,KPVT,INFO)
  37.         IF (INFO .NE. 0) GO TO ...
  38.         DO 10 J = 1, P
  39.         CALL DSPSL(AP,N,KPVT,C(1,J)) 10 CONTINUE LINPACK.  This version dated
  40.      08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab.
  41.      Subroutines and Functions BLAS DAXPY,DDOT Fortran IABS
  42.  
  43.  
  44.  
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.